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ABSTRACT 

We present a detailed study of the X-ray dust scattering halo of the black hole candidate Cygnus X-1 based 
on two Chandra HETGS observations. Using 18 different dust models, including one modified by us (dubbed 
XLNW), we probe the interstellar medium between us and this source. A consistent description of the cloud 
properties along the line of sight that describes at the same time the halo radial profile, the halo lightcurves, and 
the column density from source spectroscopy is best achieved with a small subset of these models. Combining 
the studies of the halo radial profile and the halo lightcurves, we favor a geometric distance to Cygnus X-1 of 
d = 1.81 ±0.09kpc. Our study also shows that there is a dense cloud, which contributes ^50% of the dust 
grains along the line of sight to Cygnus X-1 , located at ^ 1 .6 kpc from us. The remainder of the dust along the 
line of sight is close to the black hole binary. 

Subject headings: dust — scattering — distance — X-rays: ISM — sources: Cygnus X-1 



1. EVITRODUCTION 

It has long been known that X-ray sources with high col- 
umn densities are associated with dust "scattering halos", 
which are formed by the scattering of the sources's X-rays 
from a foreground dust containin g cloud in the in terstellar 
medium (ISM). Fi rst predicted bviOverbec^ (119651) and then 
expanded upon by lTriimper & Schonfelded (Il973h . sca ttering 
halos were first observationally confirmed by lRolj (1198 3) us- 
ing an Einstein X-ray Observatory observation of the X-ray 
binary GX 339-4. Such halos wer e later studied with virtu- 
ally all imaging X-ray i nstruments (fPredehl & Schmittll 19951: 
Predehl&Klosgfigga iPredehl et al.l 120001: iCostantini et alj 



20051) . 

The properties of a dust scattering halo, i.e., the halo ra- 
dial profile and the delay between source lightcurves and 
halo lightcurves, depend upon the composition, size distri- 
bution, and spatial distribution of the intervening, scattering 
dust grains. Therefore scattering halo s have been widely 
used to probe the porosity of g rains dMathis et al.1 [19951: 
[Smith & Dwek 1998: S mith et al.l 12002'). grain composition 
ICpstantini et al. 2005), and the spatial distri bution of the 
dust along the l ine-of-sight (LOS: IXiang et alJ [20051 \20Wi 
[Ling et alj|200^ . Comprehensive studie s of dust halos in- 
clude the systematic ROSAT studies of IPredehl & Schmitj 
( 119951) — yielding relationships between the equivalent hy- 
drogen column, A^H, and the dust grain size distribution in 
the interstellar medi um (ISM) — and the Chandra ACIS- 
S/HETGS studies of 'Xiang et al] ([2005), for the determina- 
tion of dust spatial distributions. See also lXiang et al.l ( l2007h 
where non-uniform distributions were explored in the context 
of the structure of the Milky Way. 

We note, however, that many of these results also depend 
upon the properties of the dust assumed in the modeling. 
Comparing reddening-based hydrogen column densities A^h 

[jxiang@cfa.harvard.edu; jclee@cfa.harvard.edu | 



with dust models i n the case of 4U 1724-30 7 and X Persei, 
iValencic & SmithI (120081) and iValencic et al] (2009) showed 
that o nly some of the models summarized by iZubko et alj 
(12004 hereafter ZDA) are consistent wit h the measur e ments , 
while the commonly-used models of Mathis et al.l (119771 
hereafter MRN) and Weingartner& Draine (2001', hereafter 
WDOl) significantly underestimated the measured A^h- 

The Chandra Advanced CCD Imaging Spectrometer 
(ACIS: [Garmire et alj 120031) . with its high angular and en- 
ergy resolution, is perhaps the best instrument to observe 
spatially resolved X-ray dust halos. Furthermore, via tim- 
ing analysis of the halo these data can be used to determine 
the distance to a bright, variable source . One such distance 
measurement was presented by Prede hl et alj (120001) for an 
observation of Cygnus X-3, using the method proposed by 
lTrumper"&" Schonfelder (1973). Since then, this method has 
been used to determine the distances to several X-ray sources, 
including 4U 1624- 490 (Xiang etal. 2007) and Cen X-3 
((Thompson & Rothsc hild 2009).' In the case of 4U 1624-490, 
the distance estimate to the object was improved from about 
30% uncertainty to 15%. 

Cygnus X-1 is a well known high-mass X-ray binary 
(HMXB) that includes a blue supergiant star and a black 
hole candidate. Although it was discovered in 1964, its dis- 
tance has been uncertain. Distance estimates ranged from 
2.5±0.4 kpc (derived from optical extinction measurements: 
iMargon et al.lll973[) to 14^0 4 kpc (based on multiple-epoch 
phase referenced VLBI observations: Lestrade et al. 1999). 
Most workers seem to have resorted to 2.0 ± 0. 1 kpc, which, 
e.g., has been used to estimate the position of the scattering 
dust dLing et al.ll2009l) . Very recent estimates, based upon 
parallax measurements of the radio jet, place Cygnus X-1 at 
1.86±0.12kpc (Reid et al., 2011, submitted). In Section[5] 
we present distance determinations independent of all of these 
estimates, and demonstrate that the halo method yields a dis- 
tance estimate of comparable quality to the parallax method. 
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Fig. 1 . — Geometry of the X-ray-scattering process for single scattering. 

In this paper, we first give a brief theoretical descrip- 
tion of dust halos in Section |2] Then, using two Chan- 
dra ACIS-S/Hi gh Energy Tra nsmission Gratings Spectrom- 
eter (HETGS; Canizares et al.li2005) observations (Section [3] 
SectionUJi, we present a study of the composition, density, and 
spatial distribution of the dust grains along the line of sight 
(LOS) based upon the halo radial profile (SectionlDi. We then 
give a detailed study for the distance determination based on 
halo timing analysis and show that this method is very reliable 
(Section |5]l. We summarize our conclusions in Section |6] 

2. THEORETICAL AND HISTORICAL BACKGROUND 

The theory governing the observed halo surface bright- 
ness has been well explored by many aut hors (e.g., 
Overbeck|[l9 65; Mauche & Gorenstein 1984; M athis & Leel 
19911; iPredehl & KlosQ il996i; .Smith et al.i ,2002i. and refer- 



ences therein) and the calculation of the time delay of a scat- 
tered photon with resp ect to an unsc attered one also has been 
discussed extensively (iTriimper & Schonf elder 1973). In the 
following, we give a brief summary of the main points of this 
theory required for our analyses. Readers desiring a more 
rigorous explana tion are referred to the relevant papers cited 
above. See also IXiang etaTI (120071) for details (including a 
flow chart) of the analysis sequence. 

2.1. Halo Surface Brightness 

X-ray scattering halos are formed when X-rays emitted by 
a background source a re scattered by dust i n the intervening 
ISM. A s discussed by i Mafliis & Leg (119911) and ISmith et all 
(l200l,for a single scattering (Fig.[T]l at an observed angle 9, 
the observed first-order scattering halo surface brightness for 
photons with energy E can be described by 



lili{e,E) = Fx(E)NH nia) f(x)(\-xr^ 

xS(a,E,9sca) dadx , (1) 

where Fx(E) is the total observed X-ray photon flux (in units 
of photons cm~^ s~') at energy E, the relative distance x = 
d/D IS, the ratio of the distance from scattering grain to the 
observer id) and the distance from the source to the observer 
(D), A^H is the equivalent hydrogen column density between 
the observer and the X-ray source, and f{x) is the ratio of the 
local hydrogen density at xD to the average hydrogen density 
along the LOS (Fig. [U. The grain properties are defined via 
the size distribution of the dust grains, n(a), and the energy- 
dependent differential cross section at a scattering angle of 
9sai for a spherical particle of radius a, S(a,E , 9^^:^), which in 
the Rayleigh-Gans approximation is given by (iMathis & Led 
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where ci = 1.1 x lO-'-cm^/sr, K = 0A575(E/kevf(a/^imf, 
Z is the mean atomic charge, M is the molecular weight (mea- 
sured in atomic mass units), p is the mass density and F{E) 
is the atomic scattering factor (taken, e.g., from Henke 1981). 
Thus a photon of wavelength A will be scattered into a typical 
angle 

61 = — «£"'. (3) 
TT a 

F or small angles (sever al arcmin), 9 ^ {I -x)0sca- 

iMathis & Led (119911) have shown that for scattering optical 
depths greater than 1, doubly scattered radiation may domi- 
nate the multiple scattering at distances of several arcmin. For 
the scattering halo discussed here, however, even at energies 
as low as 1 keV the optical depth, r < 1. Multiple scattering 
therefore can be neglected for the photon energies considered 
in this work. 

Eq. [T] shows that the halo surface brightness, as well the 
time delay of the scattering photon (see 32.21 below), tightly 
depends upon the spatial distribution of the dust grains and 
on the size distribution of the grains. As an example. Fig. |2] 
shows halo radial profiles for different spatial distributions. 
Halos due to dust clouds close to a source are sharp and nar- 
row, while halos due to clouds closer to the observer are flat 
and wide. 

2.2. Time Delay of the Scattering Photon 

From the geometry shown in Fig.[T] the total distance trav- 
eled by a photon scattered at xD is 

xD 



D 



sm{9) 



+ y^{D-xDy + {xD tan 61)2 



(4) 



This photon will travel a longer distance than the unscattered 
one. For small angles, 6* ^ 1 and sin(6') ^ tan(0) ~ 9. In this 
case, the delay time is 



dt = 



D -D 



1.21 sec 



D 



Ikpc J VI" 



1- 



(5) 



The delay time therefore increases dramatically at large an- 
gles and easily becomes greater than the time scales of the 
intrinsic intensity variations. For this reason, if one wants 
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Fig. 2. — Halo radial profiles for the different spatial distributi ons o f the 
dust grains based on the dust model of ZDA BARE-GR-S. (See j|2.3l for a 
discussion of this model.) The hydrogen column densities are based on the 
halo fits discussed in ^ and the indicated ranges are the relative distance 
between observer, at 0, and source, at 1 . The specific model shown in fact 
corresponds to the first entry in Table|3] 

to determine the distance from time delay measurements by 
inverting Eq. |5] the variability of the halo at small scattering 
angles has to be used. In addition, as shown in Fig.|5]the delay 
time tightly depends on the position of the scattering medium. 
For this reason, the distance and the spatial distribution of the 
dust have to be determined together 

2.3. Dust Grain Models 

As discussed above, dust composition and grain size distri- 
bution affect the overall determination of the halo properties. 
In this section, we will briefly describe the most commonly 
used models for dust in the ISM. These include the two most 
commonly used grain models of MRN and WDOl, as well 
as a sequence of 15 differe nt kinds of grain models recently 
constructed by Zubko et al. (2004). 

The classical grain model is that of MRN, which assumes 
both graphite and silicate grains with a power law distribu- 
tion. Based on the observed interstellar extinction over the 
wavelength range 0.11 /xm < A < 1 /im, the size distribution 
is described by n(a) = K{a/\ /im)"-'-^, where K = lO"''''^'* for 
graphite and K= 10"'^^' for silicate. Both size distributions 
ar e sharply cut off outside of 0. 005 ^m < a < 0.25 /im. 

IWeingartner & Draind (120011) constructed grain models for 
different regions of the Milky Way, LMC, and SMC. In con- 
trast to the MRN model, the WDOl model includes suffi- 
ciently small carbonaceous grains (including polycyclic aro- 
matic hydrocarbon molecules, PAHs) to account for the ob- 
served infrared and microwave emission from the diffuse 
ISM. The size distributions of both carbonaceous and silicate 
grains are not simple power laws (see Fig.©. 

Unlike MRN and WDOl models, which based elemen- 
tal abundan ce constraint s upon only one interstellar medium 
abundance, IZubko et"an ( |2004 ) include solar (S), F- and G- 
star (FG), and B-star (B) abundances when they derive their 
interstellar dust models. The ZDA model is derived from si- 
multaneous fits to the far-ultraviolet to near-infrared extinc- 
tion and assumes five different dust constituents: (1) PAHs; 
(2) graphite; (3) amorphous carbon; (4) silicates in the form 
of olivines (MgFeSi04); and (5) composite particles contain- 
ing different proportions of silicates, organic refractory mate- 
rial (C25H25O5N), water ice (H2O), and voids. ZDA consider 
two groups of models, one group that include only PAHs and 
bare grains and another group containing further composite 
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Fig. 3. — Example grain-size distributions for the MRN, WDOl, and the 
ZDA BARE-GR-S dust models. 

particles. These groups are called BARE and COMP, respec- 
tively. The BARE and COMP models are further subdivided 
according to the properties of carbon in the different models. 
ZDA distinguish between graphite (GR), amorphous carbon 
(AC), and no carbon (NC). Taking into account the different 
abundances (designated previously as -S, -FG and -B), a to- 
tal of 15 different dust grain models are obtained. The des- 
ignation of each model derives from the abbreviations listed 
above. For example, the model consisting of bare silicate, 
PAHs and graphite and derived by assuming F- and G-star 
abundances will be called BARE-GR-FG. 

In order to compare these three kinds of models with each 
other. Fig. [3] shows their size distributions. In order to avoid 
unnecessary clutter in the figure, for the ZDA models we only 
plot the size distribution of BARE-GR-S. This latter model 
yields the proper redd enings along the LO S to the low-mass 
binary 4U 1724-307 dValencic et al.ll2009h . We also Hst the 
key parameters, e.g., composition, size range and abundance 
of dust grain, in Table [T] The parameters of models 1 to 15 
are ta ken from (Zubko et al. 2004), MRN from (Mathis "etaLl 
Il977h . and WDOl from ( Weingartner & Draine 2001) respec- 
tively. The model of "XLNW" is a modified form of ZDA 
BARE-GR-S (listed as model 1 in the tables; see ^for a de- 
tailed discussion). 

Based on a detailed analysis of observed XAFS (X-ray Ab- 
sorption Fine Structure) in Cygnus X-1 spectra of the Fe L 
band (Lee et al. 2009a; Lee et al. 2011) using the techniques 
described in Lee et al. (2009b, see also Lee 2010), we intro- 
duce an additional model, which we dub XLNW, where the 
major Fe-containing grain compound olivine (FeMgSi04) in 
the ZDA BARE-GR-S model is replaced with an iron grain 
consisting of a metallic iron core surrounded by troilite (FeS), 
and enstatite (MgSiOa). In generating this new model, for 
each compound, the normalization coefficient (A) is adjusted 
to meet the criteria of the following dust mass equations (in 
different representations on both the left and right side) 



AzM 
Na 



dan(a)a 



(6) 



where for the grains replacing olivine, Az = A^g = Asi = Ape 
is the abundance of the dust compound, M is the molecular 
weight in atomic mass units (AMU), and A^a = 6.02 x lO^-' 
is Avogadro's constant. As described above, 
and n(a) are, respectively, mass density, radius, minimal ra- 
dius, maximal radius and size distribution of the dust grains. 
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TABLE 1 
Summary of dust models. 



Model" 


Model* 


PAH' 


Graphite'' 


ACH2'' 


Olivine^ 


Composite^ 


Iron'' 


Reference' 


No 


Name 


Size Range 


Size Range 


Size Range 


Size Range 


Size Range 










Abundance 


Abundance 


Abundance 


Abundance 


Abundance 


Abundance 




1 


BARh-CjR-S> 


0.35-5.0 


0.35-330 




0.35-370 












33.0 


212.9 




33.3 




33.3 




2 


"TJ A T~>T7 /""Tl TT/"* 


0.35-5.0 


0.35-300 




0.35-340 












35.2 


212.3 




33.1 




33.1 




3 


BAKb-CjR-B 


0.35—3.5 


0.35—320 




A T C TTA 

0.35—320 












33.3 


221.1 




28.5 




28.5 




4 


CuMP-UR-S 


0.35-5.5 


0.35-500 




0.35—440 


A AAA 

20-900 










33.5 


109.2 




25.0 


8.0 


33.0 




5 


CUMF-CjK-rCj 


0.35—5.5 


A 1 C T AA 

0.35—390 




A 1 C T AA 

0.35—390 


20-750 










35.8 


133.3 




26.1 


6.3 


32.4 




6 


CUMF-LrK-B 


0.35—5.5 


A T C CIA 

0.35-520 




0.35—330 


20-^50 










33.7 


133.0 




20.0 


7.8 


27.8 




7 


BAKb-AC-S 


0.35—3.7 




20-260 


3.5—370 












51.4 




213.6 


33.5 




33.5 




8 


BAKb-AC-rCj 


0.35—3.6 




OA O OA 

20-2 80 


3.5-370 












52.4 




212.7 


34.2 




34.2 




9 


BAKb-AC-B 


0.35—3.7 




20-250 


3.5—330 












52.2 




223.1 


28.7 




28.7 




10 


CUMF-AL-S 


0.35—3.8 




20-250 


0.35—400 


20-9 1 










50.6 




75.2 


23.7 


9.6 


33.5 




1 1 


L. U JVIF- AL.- r Lr 


0.35—3.5 




20—250 


(J. JJ— 4(JU 


zt)— odU 










51.7 




81.2 


24.5 


9.2 


33.7 




12 




yj.jj—j.y 




22—2 1 


fi ?sn 

U. J J — ZJU 


20—700 










51.5 




28.1 


14.3 


13.6 


27.9 




ID 




U. J J— J.D 






fj ac "I /in 
yJ.jj—J'+V 


ZlJ— OUU 










so n 

JU.U 






1 S 7 


1 A 7 






14 


COMP-NC-FG 


0.35-3.5 






0.35-360 


19-850 










51.0 






19.1 


14.8 


33.9 




15 


COMP-NC-B 


0.35-3.8 






0.35-180 


20-800 










51.5 






12.7 


15.4 


28.1 




16 


MRN 




5-250 




5-250 






2 








270.0 




33.0 




33.0 




17 


WDOl 




0.35-1000 




0.35^00 






3 








231.0 




36.3 




36.3 








PAHs 


Graphite 


ACH2 


Enstatite 


Fe metal 






18 


XLNW 


0.35-5.0 


0.35-330 




0.35-370 


0.35-370 




4 






33.0 


212.3 




33.3 


33.3 


33.3 





Model number and model name 

'^^S The unit of the size range is 10""* (nm) and the abundance is the molecular number of composition per hydrogen in ppm (10^^). The abundance 
of composites is based on the abundance of olivine residing in the composite. 

^ The silicate in the ZDA, MRN and WDOl models is explicitly shown as olivine (MgFeSi04) in order to discriminate it from other silicates, e.g., 
enstatite (MgSiOj) in the XLNW model. 

^ The iron abundance equals to the sum of abundances of all iron compounds, e.g., olivine, composite and iron metal. The unit is also ppm (10^^). 
' References for dust model parameters: 1 . IZubko et"an )2004ll , 2. IMathis et alJll97TI) . 3. IWeingartner & Drai"n3 120011) and 4. this paper 



(The normalization coefficient. A, is a parameter contained 
within n(fl).) For the iron compound, i.e., iron metal + troilite 
(FeSo.e), p = 5.2 gcm-^ Az = 33.3 x 10"'' and M = 75.2 AMU 
while p = 3.2 gcm-^Az = 33.3 x lO-*" andM= 100 AMU for 
enstatite (MgSiOs). Values for p are taken from a mineralogy 
databaseQ; where the compound is an admixture, p is taken to 
be the weighted average of the minerals which make up the 
compound. 

As in the ADA and MRN models, all Fe and Si are assumed 
to reside in the dust, as borne out also by our XAFS fitting. 
The particle size distributions of iron dust and enstatite are 
the same as those of silicate in ZDA BARE-GR-S. The nor- 
malization of the size distribution, however, has been slightly 
changed by adjusting the amount of dust appropriate for the 
assumed depletion. The size distributions and the normaliza- 
tion coefficient (A) of the XLNW dust model are shown in 
Fig. El 



http://webmineral.com/ 



3. OBSERVATION AND DATA ANALYSIS 

We now will use two Chandra -HETGS observations to 
study the X-ray halo of Cygnus X-1 using the methods dis- 
cussed above. The first observation occurred on 1999 October 
19 (ObsID 107) and lasted about 15 ks. During this time the 
source was transiting between its hard, po wer-law dominated 
state, and the more thermal soft state (see iSchulz et al.ll2002l 
for a study of the source behavior). As is characteristic for 
black holes in this transitional state, the overall source vari- 
ability was low. This property allows us to use these data to 
measure a single non-variable radial halo profile which we use 
to determine the spatial distribution of the dust grains. The 
second observation was performed on 2004 April 19 (ObsID 
3814). This later observation lasted about 50 ks and was de- 
signed by us to be performed during the upper conjunction of 
the black hole (Hanke et al. 2009). During this orbital phase 
the X-ray lightcurve is strongly variable due to frequent pho- 
toabsorption dips as the line of sight to the black hole crosses 
through the clumpy focused stellar wind of the companion. 
Due to dramatic variability, this observation is well-suited to 
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Fig. 4. — The particle size distributions of our XLNW dust model, which FiG. 5. — The streak spectrum and MEG first order spectra overplotted with 

is a modified form of the ZDA BARE-GR-S model. the best fit model. 



determine the distance to Cygnus X- 1 . 

We create two datasets, as described below, to perform 
these analyses. One is a steady spatial/spectral halo radial pro- 
file that is used for determining the composition of the halo. 
This dataset is primarily comprised of Obsid 107, however, 
we describe a procedure to replace the inner (<20 arcsec) data 
in order to obtain a more accurate profile. The second dataset 
consists of variable halo lightcurves, as a function of radius, 
comprised solely of data from ObsID 3814. The latter dataset 
is used for the distance determination. 

In order t o form the da tasets, we follow the same procedure 
described in lXiang et alj (j2007) to extract the HETGS spectra, 
halo radial profile, and halo lightcurve. (See that paper for a 
detailed analysis flow chart.) CIAO 4.2 with CALDB 4.2 was 
used. A theoretical point spread function (PSF) created with 
the Chandra Ray Tracer (ChaRT; Carter et al. 2003) is used 
for small angles (<60 arcsec). Cha RT has been shown to rep- 
resent properly the PSF behavior; ( IXiang et al.ll2007b . The 
PSF obtained from ChaRT simulations, however, is underes- 
timated at large off-axis angles (> 60"). We therefore use the 
bright and halo-less point source Her X-1 to perform an ab 
initio determination of the PSF at large angles. Her X-1 was 
observed on 2002 July 02 (ObsID 3662) with Chandra ACIS- 
I. The parameters of these three observations are listed in Ta- 
ble|2] 

Both the halo and the PSF represent functions of radius and 
energy. The PSF specifically is the ratio between the sur- 
face brightness, absent a halo, and the point source spectrum. 
In order to determine the surface brightness associated with 
the combined halo and PSF, we need an accurate estimate of 
the source spectrum for both Cyg X-1 observations discussed 
here. In each of these cases, the spectra extracted from the 
zeroth-order suffer heavy pileup over the entire energy band 
and also suffer pileup in the first order gratings s pectra above 
0.7 ke V. We therefore follow the method used by lSmith et alJ 
(120021) and lValencic et al.l (120091) to extract the spectrum from 
the readout streak, which is completely pileup free. We then 
use the simple model - a disk blackbody plus powerlaw cou- 
pled with cold gas absorption - that has been used to fit these 
two observations ( Schulz et al. 2002; Hanke et al. 2009). Be- 
low 0.7 keV, we use the MEG±1 spectra, which are pileup 
free. The flux of Cygnus X-1 is determined from the fit 
to these joint MEG±1 spectra and the readout streak spec- 
tra. The transitional state spectra, overplotted with the best 
fit model, are shown in Fig. |5j (The spectral flux during the 



hard state is also derived with the same methods.) It should 
be noted that the A^h fit from these spectra contains gas local 
to the Cygnus X-1 system and therefore is expected to result 
in A^H values higher than those that we expect to obtain from 
modeling the halo radial profile. 

As discussed above, ObsID 107 was used to generate the 
non-variable halo radial profiles. Due to the brightness of 
Cygnus X-1 during the transitional observation (ObsID 107), 
however, these data suffer severe pileup at small angles. In or- 
der to investigate clouds very near the source, halos at small 
angles are needed. We therefore excise the core of the halo 
from ObsID 107 and substitute data from the hard state ob- 
servation, where due to its low flux, pileup free halos can be 
extended to as low as 7" radius. In this work, the total sur- 
face brightness (see Fig.|6l) at greater than 20 arcsec radius is 
extracted from the transitional state observation (ObsID 107), 
while the surface brightness at less than 20 arcsec radius is 
from the hard state (ObsID 3814). Even though Cygnus X-1 
is variable, this approach is justified since at small angles the 
time delay between the source and the halo is negligible, and 
we therefore can compare the time averaged source flux to the 
time averaged radial profiles. 

To determine the energy dependent halo profile using these 
data, we extracted the halo and PSF radial profiles in energy 
bands of 200eV width over the range 0.4-10. OkeV. That is, 
we use bands that correspond roughly to the energy resolution 
of ACIS-I/S, which varies from 100 to 200 eV in this range. 
We find that above 4.0 keV, the halo intensity is much smaller 
than the PSF intensity (i.e., the surface brightness attributable 
to the point source). Additionally, the Rayleigh-Gans ap- 

J iroximation underesti mates the halo intensity below l.OkeV 
Smith &Dwe3ll998l) . We therefore only use the halo from 
l.OkeV to 4.0 keV in the following. The total surface bright- 
ness (Halo -I- PSF), PSF, and PSF-subtracted net halo profiles 
are shown in Fig. |6] for the 2-2.2 keV band. In this band, at 
angular distances greater than 10", the halo surface bright- 
ness is significantly larger than the surface brightness of the 
point source (i.e., the PSF intensity). On the other hand, for 
all energies >1.6keV and angular distances < 10" the data 
are fully dominated by the point source and therefore are ig- 
nored in the following study. For energies below 1.6 keV we 
can approach the point source to separations as low as 7". At 
smaller angular separations the pileup ratio is greater than 5% 
and the data have to be ignored at all energies. 
A final complication arises for energies above 2.0 ke V and 
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TABLE 2 

Parameters of three Chandra observations 



Source 


Observation 


State 


Instrument 


Start Time 


End Time 


Exposure 


Usage" 




ID 










(ks) 




Cyg X-1 


107 


Transition 


AICS-S/HETG 


1999-10-19 19:17:37 


1999-10-19 23:52:10 


11.4 


RP 


Cyg X-1 


3814 


Low/Hard 


AICS-S/HETG 


2003-04-19 16:47:31 


2003-04-20 06:41:42 


47.2 


RP + LC 


Her X-1 


3663 


High/Soft 


AICS-S/I 


2002-07-01 23:37:38 


2002-07-02 14:03:05 


49.6 


PSF 



^ "RP" means that the data are used to generate the^Halo Radial Profile". "LC" 
respectively. 



for "Halo Lightcurves" and "PSF" for^Point Spread Function", 



Total Surface Brightness 
PSF Surface Brightness 
Halo Surface Brightness 




50 

Radius (arcsec) 

Fig. 6. — Comparison of the radial profiles between the halo and the PSF 
(i.e., the surface brightness attributable solely to the point source) in the 2.0- 
2.2 keV energy band. The total surface brightness (black) was created from 
a combination of ObslD 3814 «20 arcsec) and ObsID 107 (>20 arcsec). - 
A vertical black dashed line delineates these two regions. The PSF profile 
was created from ChaRT simulations (<60 arcsec) and from a non-variable 
Her X-1 observation (> 60 arcsec). Data from less than 10 arcsec (marked 
by vertical blue dashed line) is not considered in this energy band. 

radial distances > 200". Here, the halo profile is contami- 
nated by photons from the P' order gratings spectrum. Of 
special importance are photons in the MEG±P' order, which 
starts closer to the zeroth order than the HEG±P' order spec- 
trum. We therefore ignored the halo radial profile at > 200" 
for energies above 2.0 keV. 

4. SPATIAL DISTRIBUTION OF DUST GRAINS 

For our initial modeling, dust uniformly distributed along 
the entire LOS was assumed. This approach did not result 
in satisfactory descriptions of the halo radial profile observa- 
ti ons using any of the 1 8 du st models. Following the method 
of lXiang et all(l200ll2007l ). we therefore assume an inhomo- 
geneous dust distribution along the line of sight. We divide 
the LOS into 10 approximately logarithmically spaced bins. 
That is, the spatial resolution of the dust distribution increases 
towards the source. The relative density distribution along the 
LOS as derived from fits to the halo radial profile and assum- 
ing the ZDA BARE-GR-S model is shown in Fig. [T] The 
ZDA BARE-GR-S model gave the best fits for the reddening 
of 4U 1724-307 in the work by Valencic et al. (2009) and it 
is also one of our favored models (see Section|6]l. 

The data prefer approximately three regions of different 
characteristic densities. The highest density region is found at 
a relative distance between « 0.84-0.94. The regions extend- 
ing from K, 0.94-1.0 and from w 0.0 — 0.6 exhibited lower 
columns. Similar results were obtained when using 10 lin- 
early spaced bins (inset of Fig. |7]i: there was little column 
between 0.0-0.8, and the bins between 0.8-0.9 and 0.9-1.0 
have different characteristic densities. 

In the ensuing analysis, we discuss these three characteris- 
tic regions as if they were three individual clouds. Further- 
more, we label these three regions, starting near us and work- 



ed 



0.9 0.8 
Relative Distance x 



0.6 



0.0 



Fig. 7. — Relative density distributions along the LOS based on halo radial 
profile fits using the ZDA BARE-GR-S dust model. The observer is at 0. and 
the source is at 1. A cloud at ~0.88 is obvious. The inset (a) figure (in blue) 
also shows a cloud near the source when ten linear-equally spaced bins are 
instead used in the model. (In subsequent discussion and fits presented below, 
a simplified three-zone model, inspired by the above results, is used.) 

ing outwards towards the source, as Cloud Obs, Cloud Mid, 
and Cloud Src. Cloud Mid contributes more than 50% of the 
total dust g rains along the LOS. This result is consistent with 
iLing et al.l (12009). who use a cross-correlation method to de- 
rive the position of this dominant cloud from the delay time 
of the scattered photons. 

Figure |2] shows that the halo flux between 20" and 45" 
is dominated by photons scattered by Cloud Mid. The de- 
lay time of scattered photons depends especially strongly on 
the distance of this cloud. To make the calculation of the 
predicted distances and fitting of the halo lightcurves more 
tractable, we further simplify the description of the dust spa- 
tial distribution. Specifically, based upon the preliminary ra- 
dial profile fit shown in Fig. |7] that exhibits three regions 
with different characteristic densities, we formally divide the 
LOS into three components. Uniform density regions are 
used to represent Cloud Obs between 0.0- Xn, Cloud Src be- 
tween x,- 1 . 0, and Cloud Mid centered at Xc . iLing et aP (120091) 
claimed limits on the width of this latter, dominant cloud of 
Ax = 0.016;!;Q 0Qg; therefore, in this work we have adopted a 
fixed width of 0.016 (^30 pc). We reexamine this limit in 
^ where we show that a cloud with a large extent does not 
fit the dust scattering halo lightcurves. For convenience, we 
take the upper limit tohe x= 1 ; however, it should be noted 
that the dust located very near the source, i.e, x > 0.99, might 
be contaminated by the wind of HDE 226868, the companion 
of Cygnus X-1 . Note that the halo scattered by this dust is 
concentrated at very small angles (<7") and contributes very 
little to the portion of the halo considered in this work. We 
therefore still use the same dust model along the entire LOS. 
We have also checked that changing the upper limit of Cloud 
Src to X = 0.99 from x = 1 .0 does not provide significantly dif- 
ferent results for the A^h due to the Cloud Src nor does it affect 
the distance determination. 
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Fig. 8. — Point source lightcurves and lightcurves at a given annulus, in 
tlie l^keV band, overplotted by tlie best-fit models (red) using the ZDA 
BARE-GR-S dust model. The blue line indicates the lightcurves derived from 
the dust model, while the green line is the modeled PSF lightcurve. Their 
sum yields the red lightcurve. The top panel (a) is the central point source 
lightcurve. The middle panel (b) is the 20-25" annulus. The bottom panel 
(c) is the 30-35" annulus. 

The relative distances Xq, Xq, and were initially allowed to 
vary in our fits. We found that both the total hydrogen column 
density and the position of Cloud Mid, i.e., Xc, are not greatly 
affected by variations of Xq but are strongly correlated with 
the value of Xs. The former fact is not surprising given the low 
relative density in this region close to the observer, as shown 
in Fig. |7] For the following halo radial profile fits we there- 
fore fix Xo to 0.70, which is the best fit value for the model 
of ZDA BARE-GR-S. On the other hand, x^ was allowed to 
vary between 0.7 5-0.9 8. The eighteen dust grain models de- 
scribed in Section O— MRN, WDOl, 15 ZDA models, and 
our XLNW model — are used to fit the halo radial profile. 

When the XLNW dust model is used to fit the halo radial 
profiles, we derive a scattering hydrogen column density that 
is consistent with the absorption hydrogen column density de- 
rived from XAFS studies (Lee et al. 201 1). It should be noted 
that our XLNW model yields a similar position for Cloud Mid 
as do the BARE-GR-S, and BARE-GR-FG. These latter mod- 
els also yield scattering hydrogen column densities consistent 
with the XAFS studies. 

We now turn to using the above preliminary fits to the halo 
radial profile to constrain the fitting of the halo lightcurves. 

5. DISTANCE DETERMINATION 

The halo lightcurves of Cygnus X-1 clearly show signifi- 
cant delays with respect to the source lightcurves. The delay 
time strongly depends on the observed off-axis angles. In or- 
der to measure the delay, we calculate halo light curves for 
six, 5" wide annuli, with the inner radii ranging from 20^5". 
Outside this interval, the delay times exceed the duration of 
our observation. 

Following the method of Xiang et aTl (l2005h . we use the 
derived dust spatial distribution to fit the halo lightcurves with 



the absolute distance to the source being the free parameter. 
The best-fit distance is determined by minimizing based on 
simultaneously fitting the delays in all lightcurves. Figure [8] 
shows the point source and halo lightcurves overplotted with 
the best-fit of the ZDA BARE-GR-S model, using a single fit 
distance for all lightcurves considered simultaneously. 

The time delayed lightcurves also carry information about 
the positions of the clouds such that the lightcurve fits also 
can be used to constrain the position of x, further than by just 
using the halo profile. We therefore vary Xj from 0.75 to 0.98, 
and at each trial value of Xj we perform a constrained fit to 
the radial profile to determine x^ and the hydrogen column 
densities. We then use these radial profile derived parameters 
to fit the halo lightcurves and obtain a minimum x^- While the 
from solely a fit of the halo radial profile varies little with 
Xs, the x^ from fits of the halo lightcurve varies significantly 
with Xs. We use this variation in x^ to create a lightcurve 
determined uncertainty for Xs (i.e., Ax^ = 2.71 to determine 
the 90% error bars). For 13 of the 18 dust grain models, the 
values of Xs derived from the preliminary radial profile fits 
described above and those values derived from the lightcurves 
were self-consistent with one another, i.e., the 90% error bars 
for each estimate (lightcurve, radial profile) of Xs significantly 
overlapped. 

Given that the iterative approach to fitting the radial profile 
and then fitting the lightcurve failed to produce self-consistent 
fit parameters for five cases, we then turned to a joint fit of 
these data. Specifically, we used the solutions obtained from 
the iterative approach as starting parameter values, and then 
simultaneously fit the radial profile and the lightcurve data 
with six variable parameters: Xs, Xc, d, and three hydrogen 
column densities, A^^'"*, and A^^'^ The results of these fits 
are presented in Tabled The N"^' are the sum of A^™'', 
and N^'^. The iron column density A^p°' is derived from N^^ 
based on the iron abundance for each dust model (see Table[T|. 
Parameter uncertainties are given as the 90% confidence level 
for one interesting parameter, i.e., Ax^ = 2.71 for fits holding 
the parameter of interest frozen at a given value. Values of Xj, 
Xc, A^p°', and d for each dust model are presented in Fig.|9l 

Our derived scattering dust iron column density varies from 
0.94 X 10" cm"2 to 1.92 x 10'^ cm"^ which should be equal 
to the Galactic absorption dust iron column density. We there- 
fore use results from absorption spectra to constrain further 
the dust models. We choose to compare to the iron col- 
umn density as it can be measured directly via mode ling of 
the F e L edge in high s pectral resolutio n data ( Hanke et alj 
2009t iLee et all l2009al: Lee et alj I20TTI) . As discussed by 
Wilms et alJ (l2000h . neutral column absorption in the X-rays 
is dominated by Z > 2 elements, i.e., not hydrogen. Compar- 
ing directly to the X-ray derived Fe column therefore removes 
the ambiguity of the assumed elemental abundances relative 
to hydrogen in such spectral models. 

The total absorption iron column density of Cyg X-1 de- 
rived from absorption spectra varies from 0.96 x lO" cm~^ 
dJuettet al.ll2006 ) to 2.52 x 10'^ cm'^ (HankeeLalJ|2009'). 
Aside from possible instrumental differences, these differ- 
ent values are in part due to variations in the column den- 
sity local to the source, e.g. , because of absorption in the 
stellar wind (see iHanke et al.l l2009). as well as variations in 
the assumed absorption cross sections. Our recent studies 
of observed X-ray Absorption Fine Structur e (XAFS) in the 
Low/Hard, Transition, and High/Soft states (iLee et al.ll201 ll) 
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Fig. 9. — The results of the joint fits to the halo radial profile and lightcurve for each model. (Model ID is shown along the .t-axis, and model names are given 
in the top panel.) The top two panels give the relative positions of the cloud nearest the source (Xs), and the next closest cloud (.Vc). The next panel gives the 
summed total neutral column (converted to an equivalent Fe column based on the respective olivine abundance listed in Table 1 for each dust model) for all of the 
clouds between us and the source. The two red lines show the iron column estimate obtained from the Lee et al. 1 20 II ) modehng of the ~700 eV FeL absolution 
edges and associatd XAFS in HETGS spectra of Cygnus X-I . For all panels, models plotted in red agree with the neutral column density estimate based on iron. 
The bottom panel gives the best fit source distance for each model. The median value is 1.82 kpc, while three of the six models that agree with the neutral column 
estimates of Lee et al. ( 201 1 .) give distance estimates that lie in the range from 1. 72-1. 90 kpc. 



suggest that the dust iron column density associated with Cyg 
X-1 is 1.62 ±0.06 X 10'^ cm"^. Only six of the eighteen dust 
models considered for the halo modeling have A^pe consis- 
tent with this value. (The parameters for these models are 
plotted in red in Fig. |9]) Of these six models, three models, 
ZDA BARE-GR-S, ZDA BARE-GR-FG, and our proposed 
model XLNW, yield source distances that lie in the range 
from 1.72-1. 90 kpc (The median distance for all three mod- 
els is 1.82 kpc), whereas three models, ZDA ZDA COMP- 
NC-S, ZDA COMP-NC-FG and ZDA COMP-NC-B, yield 
larger distances of 3.15-3.59kpc. The large distance values 
for these latter three models, compared to the most recent dis- 
tance estimates (discussed below) leads us to discount these 
fits. This leaves us with three models that we favor to describe 
the halo and absorption profile of Cygnus X-1 . 
Lastly we consider the possibility that the cloud with a rel- 



ative distance between 0.8 and 0.9 is a thick cloud, or that 
the cloud might be divided into 2 parts. For the first possi- 
bility, we set a cloud with a width of 0.1 (relative distance) 
located at = 0.878. This location is the best fit value of Xc 
derived from radial profile fits using the ZDA BARE-GR-S 
model. The distances then fitted from halo lightcurves at dif- 
ferent angles are then not consistent with each other, as shown 
in the middle panel of Fig. [10] For the second possibility, we 
locate one thin cloud at 0.85 and another thin cloud at 0.90. 
The width of each cloud is set to 0.016. The distances fitted 
from individual halo lightcui-ves are again found to be incon- 
sistent with each other, as shown in the bottom panel of Fig. 
[TOl We also noticed that the x^/DoF of the distance fit with 
thick cloud (xVDoF = 1003/777) or two thin clouds (x^/DoF 
= 983/776) is much worse than that from one thin cloud fit 
(X^/DoF = 884/777). This further justifies our assumption of 
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Fig. 10. — The distances derived from lightcurves at different radii for one 
thin cloud (a), one thick cloud (b), and two thin clouds (c). Only the distances 
derived from one thin cloud are consistent with each other for each radii. The 
dashed horizontal lines indicate the average distance. The uncertainties in 
this figure are only from the fits to the halo lightcuves. The error bars in this 
figure are only Itr. Additional uncertainties from uncertainties in the radial 
profile fits are not included. 

a single thin cloud at the position x^. 

6. DISCUSSION AND SUMMARY 

Using the X-ray dust scattering halo radial profiles and 
lightcurves combined with the X-ray absorption hydrogen 
column density, we can begin to distinguish among the dust 
grain models. Our favored dust models are BARE-GR-S, 
BARE-GR-FG, and XLNW. It should be noted that the first 
two models, i.e., BARE-GR-S and BARE-GR-FG, also yield 
the proper optical e xtinction for the LO S to the X-ray bi- 
nary 4U 1724-307 ( Valencic et 10120091) . Our fits using the 
classic MRN model and the more recent model of WDOl 
underestimated the hydrogen column density to Cygnus X- 
1 ; similar underestimates w ith these two models were found 
bv IValencic & SmithI (l2008h for X Persei and Valenc ic et all 
(120091) for 4U 1724-307. All three of our favored models con- 
tain graphite, while none of the dust models containing amor- 
phous carbon simultaneously yield good neutral columns and 



source distances (see below). This result suggests that carbon 
in the ISM prefers to reside in graphite instead of amorphous 
carbon dust. Our proposed model XLNW, which separates 
iron from olivine, can yield a proper hydrogen column den- 
sity and a proper position for the source, as discussed below. 

Our results show that the dust scattering halo radial pro- 
files together with the halo lightcurve provide a potentially 
powerful tool to determine the relative position of the cloud 
along the LOS, and to determine the geometric distance to 
the point source. All three of our favored models show that 
Cloud Src is located at --0.92-0.95 and that Cloud Mid 
is located at Xc « 0.885. Our fitted distances to Cygnus X- 
1 for the most part lie within a narrow range of 1.72-L90. 
This distance range is consistent with recent radio parallax 
measurements for Cygn us X-1 whi ch yield a distance of 
d= L86± 0.12 kpc (Reid et al.ir2011l) . This is considered the 
most reliable method to determine the distance. The distance 
derived from our proposed dust model containing iron metal 
and troilite is consistent with those derived from BARE-GR- 
S and BARE-GR-FG, both of which yield the proper column 
density. 

The first estimate of the distance to Cygnus X- 1 was made 
by Margon et al. ( 1973) using the surveyed optical extinction 
in the field immediately surrounding Cygnus X-1 . The dis- 
tance of li = 2.5 ± 0.4 kpc was confirmed by iNinkov et alJ 
(1987) and was used by man y researchers. Russell et al] 
(120071) and iGallo et al.1 d2005h imply that Sharpless 2-101 
hereafter Sh 2-101, which is is a bright reflection nebula 
to the north east o f Cygnus X-1 , is close to Cygnus X-1 . 
iRussell et al.l (l2007h also suggest that the nebula interacts with 
the jet of Cygnus X-1 . The distance of Sh 2-101 is generally 
assumed to be larger than 2.2 kpc. Our results, most of which 
are significantly lower than 2.2 kpc, hint that Sh 2-101 is a 
background source and therefore the Cygnus X-1 jet is not 
interacting with it. Alternatively, measurements of the dis- 
tance to Sh 2-101 may be overestimates. 

The authors would like to acknowledge useful conversa- 
tions with Randall Smith. J.C.L gratefully acknowledge 
funding support from Chandra grant AR8-9007, and the 
Harvard University Faculty of Arts and Sciences. M.A.N 
was supported by NASA grants GO8-9036X and SV3- 
73016. J.W was partly supported by the European Commis- 
sion under contract ITN 215212 "Black Hole Universe" and 
by the Bundesministerium fiir Wirtschaft und Technologie 
through Deutsches Zentrum fiir Luft- und Raumfahrt grants 
50OR0701 and 50OR1005. 



REFERENCES 



Canizares, C. R., et al., 2005, PASP, 117, 1144 

Carter, C, Karovska, M., Jerius, D., Glotfelty, K, & Beikman, S., 2003, 

ADASS XII, ASP Conference Series Vol. 295, eds. H.E. Payne, R.I. 

JEdzrejewski, & R.N. Hook, 477. 
Costandni, E., Freyberg, M. J., & Predehl, R 2005, A&A, 444, 187 
Gallo, E., Fender, R., Kaiser, C, Russell, D., Morganti, R., Oosterloo, T., & 

Heinz, S. 2005, Nature, 436, 819 
Garmire, G. R, Bautz, M. W., Ford, R G., Nousek, J. A.. & Ricker, Jr, G. R., 

2003, in Society of Photo-Optical Instrumentation Engineers (SPIE) 

Conference Series, ed. J. E. Truemper & H. D. Tananbaum, Vol. 4851, 28 
Hanke, M., Wilms, J., Nowak, M. A., Pottschmidt, K., Schulz, N. S., & Lee, 

J. C. 2009, ApJ, 690, 330 
Henke, B. L. 1981, in American Institute of Physics Conference Series, 

Vol. 75, Low Energy X-ray Diagnostics, ed. D. T. Attwood & 

B. L. Henke, 146-155 



Juett, A. M.. Schulz. N. S., Chakrabarty. D., & Gorczyca, T. W. 2006, ApJ, 
648, 1066 

Lee. J. C, 2010, Space Science Reviews, 157, 93 

Lee, J. C, Ogle, R M., Canizares, C. R., Marshall, H. L., Schulz, N. S., 

Morales, R., Fabian, A. C, & Iwasawa, K., 2001, Astrophys. J., Lett., 

554, L13 

Lee. J. C, et al., 2009a, in Astro2010: The Astronomy and Astrophysics 

Decadal Survey, Vol. 2010, 178 
Lee, J. C, Xiang, J., Ravel, B., Kortright, J., & Flanagan, K., 2009b, ApJ, 

702, 970 

Lee, J. C. et al. 2011, in preparation. 

Lestrade, J., Preston. R. A., Jones, D. L., Phillips, R. B., Rogers, A. E. E., 

Titus, M. A., Rioja, M. J., & Gabuzda, D. C. 1999, A&A, 344, 1014 
Ling, Z., Zhang, S. N.. Xiang, J., & Tang, S. 2009, ApJ, 690, 224 
Margon, B., Bowyer, S., & Stone, R. P S. 1973, ApJ, 185, LI 13 
Mathis, J. S., Cohen, D., Finley, J. R, & Krautter, J. 1995, ApJ, 449, 320 



10 



Xiang et al. 



TABLE 3 

Distances, cloud positions, and Wh from global fits. 
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Boundary of the region that extends from to 1 , where is the observer position and 1 is the source positon. 

^ Region midpoint. Region extends over Ac it 0.008. 
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